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1.  INTRODUCTION 

This  report  describes  the  continuation  of  the  investigation  into 
the  behaviour  of  structures  under  internal  resonance  conditions. 

Completion  of  the  analysis  on  the  three  mode  interaction  case, 
outlined  in  the  previous  report  (September  1975),  yields  the  variational 
equations  describing  the  slowly  varying  amplitudes  and  phases  of  each 
mode.  Computer  simulation  (C.S.M.P.)  of  these  equations  has  been  carried 
out.  Both  parametric  and  nonlinear  terms  are  found  to  be  of  importance 
in  the  response  of  the  structure. 

Further  experimental  work  on  the  original  four  mode  model,  where 
’parametric'  excitation  of  two  of  the  modes  leads  to  eventual  participation 
of  all  the  modes,  is  presented  along  with  a theoretical  approach  of  a more 
general  nature.  Attention  is  then  turned  to  a less  complex  two-mode 
system  exhibiting  both  parametric  and  nonlinear  behaviour. 


2.  3 MODE  INTERACTION 


Taking  the  simplified  three  mode  model  (Fig.  1),  as  described  in 
the  previous  report,  and  transforming  from  the  generalised  coordinates 
W^,  V^,  S,  A to  normal  coordinates  £ via  the  modal  matrix 
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results  in  the  equations  of  motion  in  normal  coordinates. 
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and  ap,  Xf  are  evaluated  from  system  constants. 


Upon  setting  = wocosfit,  introducing  a small  parameter  e = w0/^2> 

nondimensionalizing  with  respect  to  wq,  £ = p /wq  and  retaining  quadratic 

nonlinearities  only, 

• •• 

f(N,  N,  N,  t)  becomes:- 
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C11«2cosflt  - C12e£2(N2N*2  + M2H4  + H4N2  + N4N4  + h22  + N42  t 2N2N4) 


C2ie£2n2cosQt(N2+N4)  " C2ieVNl+H3)(N2+IV  + C22 Cl2 

• • •• 

(^l^i  + a3f,3)(a2^2  + aq^4 ) + iUtejNi  + a3N'3)(a2N2  + a^)  + 
(a2N2  + a4N4)(a11,1  + a3N3)))  " C23cZ2^2  + N4^alNl  + a3N3^ 


C^cosflt  + C32£l2(H2  + N4)(a2N2  + a^)  - C33e£2((N2  + N4) 
(N2  + N4)  + (N2  + N4)2) 


C41eA2*N2  + N4^aLHl  + a3N3^ 


Where  C^j  are  system  constants 


The  equations  of  motion  are  now  in  normal  mode  form.  Struble's  Asymptotic 
method  [1]  is  applied  in  order  to  obtain  an  approximate  solution. 

Taking  solutions  of  the  form 

Nf  = A^  (t)cos(uK  t + <^(t))  + ea^(t)  + e2b^(t)  + ... 


imposing  the  internal  and  external  resonance  conditions 


Wl  + U)2  = o>4  and  = tu4 

and  limiting  the  analysis  to  one  of  the  first  order  in  E results  in  the 
so-called  'variational  equations'.  These  are  simultaneous  first  order 
nonlinear  differential  equations  describing  the  slowly  varying  amplitudes 
and  phases  (A^  and  4>^)‘ 

In  this  case  they  take  the  form:- 
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"2Al“l 

-2A2cj24<2 


-«A„u> 


2 2 


-2A4^4>4 


-2A4W4 


C^cos^  - C2A2A4cos(<}>2  + 4>u  - 4>1 ) 

C^sin^  - C4A2A4sin(<})2  + <J>4  " 4^)  + 

C5^2A4C°S(^4  + ^ ~ C6AlA4COS^2  + ^4  ~ ^1^ 

C7fi?A4sin(<|)4  + 4>2)  - cgA1A4sin((J)2  + 4>4  - 4^)  + 2eC2u>22A2 

cgf22 A2cos (4>4  + 4*2 ) " C10A1A2COS^2  + ^4  " ^1^ 
Clin2A2sin(<^4  + <J>2)  ” Ci2Y 2sin(<t>2  + ^4  " + 


where  are  constants  for  a given  e,and  £n  are  introduced  modal  damping 
coefficients. 

* The  terms  involving  phase  angle  4>j  only  arise  from  direct  forcing 

terms  in  the  equations  of  motion,  whilst  the  terms  involving  phase  angles 

(4*0  + 4>,  “ 4>,  ) and  (<J>.  + <J>0 ) arise  from  nonlinear  and  parametric  terms 
Z 4 J.  M Z 

respectively. 

The  standard  procedure  now  is  to  set  the  left  hand  sides  of  the 
variational  equations  to  zero  to  give  steady  state  response  amplitudes  [2]. 
However  an  analytical  solution  to  the  equations  appears  to  be  unobtainable 
and  so  recourse  is  made  to  their  numerical  integration. 

On  pages  5 , 6 are  the  results  of  integration  of  both  the  complete 
equations  of  motion  and  their  corresponding  variational  equations  for  the 
model  with  exact  internal  and  external  frequency  relationships.  (cj4  = 6.20  Hz, 
u3  = 6.93  Hz,  (u2  = 22.26  Hz,  0)^  = ft  = 28.46  Hz).  Damping  coefficients, 
obtained  from  logarithmic  decrement  tests  on  the  experimental  model,  have 
been  introduced  through  a Rayleigh  Dissipation  Function.  The  excitation 
amplitude  is  0.1  mm. 

Although  the  transient  behaviour  is  different,  the  variational  equations 


do  show  the  mode  interaction  taking  place  and  the  convergence  to  steady 
state  conditions  of  comparable  amplitude  in  each  mode. 


CSMP  SIMULATION  quadratic  eauations  with  dcmi 
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3.  4 MODE  INTERACTION 

The  equations  of  motion  for  the  four  decree  of  freedom  model  (page  15 
of  previous  interim  report)  contain  both  nonlinear  and  parametric  terms, 
and  in  general  these  will  always  coexist  in  a structure.  Indeed,  the 
conditions  which  make  the  nonlinear  terms  important  in  the  behaviour  of 
the  structure  are  also  the  conditions  under  which  the  parametric  terms 
become  important,  as  shown  in  the  previous  section. 

Experimental  Work 

Further  experimental  work  on  the  four  mode  model  has  shown  that  it 
is  possible  to  excite  the  two  symmetric  modes,  fuselage  bending  'P^) 
symmetric  tail  bending  (P2)  by  forcing  'parametrically'  at  the  combination 
frequency  ft  = Growth  of  the  fourth  mode  then  leads  to  interaction 

with  the  first  and  third  modes  through  the  internal  resonance  condition 
0)((  = Wg  + <jJ^.  Hence  we  have  an  interesting  situation  in  which  the  structure 
is  excited  harmonically  far  away  from  any  natural  frequency  but  responds 
with  large  amplitudes  in  all  its  modes. 

Signal  conditioning  equipment  has  been  custom  built  to  show  the 
behaviour  of  the  model  more  clearly.  The  model  is  monitored  by  an 
accelerometer  on  the  fuselage  under  the  fin  along  with  the  strain  gauge 
pairs  at  the  roots  of  each  beam  forming  the  ' T ' tail.  The  signals  from 
these  represent  the  generalised  coordinates  of  the  mathematical  model 
used  to  generate  the  equations  of  motion.  Symmetric  and  antisymmetric 
tailplane  motions  are  monitored  by  summing  and  differencing  the  two  tail- 
plane  signals  using  analogue  devices. 

It  was  found  that  the  accelerometer  signal  and  the  sum  tailplane 
signal  were  fairly  representative  of  the  two  symmetric  modes  and  P1(. 
However,  the  fin  and  antisymmetric  tailplane  signals  required  further 
manipulation  to  give  signals  representative  of  the  antisymmetric  modes 
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and  P^.  Effectively  the  analogue  devices  do  the  same  transformation 
as  the  modal  matrix.  Fig.  2 shows  the  schematic  layout  of  the  equipment. 


The  variable  resistors  are  adjusted  until  the  signals  P,  and  P_  contain 


the  single  frequencies  to.  and  u >0  only.  These  are  then  output  to  an  ultra 


violet  recorder 


The  results  on  page  9 show  the  U-V  record  obtained  when  the  excitation 


frequency  is  32.4  Hz  (6.9  + 25.5) 


The  growth  of  the  symmetric  modes  and  eventual  interaction  in  all  the 


modes  can  be  clearly  seen.  Quasi-steady  state  conditions  with  a beat 


period  of  4.2  seconds  are  observed 
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Theoretical  Approach 


The  simplified  mathematical  model  developed  for  the  three  mode 
interaction  problem  is  inadequate  in  this  case  as  it  lacks  the  parametric 
terms  on  the  two  symmetric  inodes.  Hence  a more  general  approach  has 
been  adopted. 

The  equation 

ii  + ft  2p  = (ft  2 - ui  2 )p  + e{k  -p.cosft.et  + 1 ..p.p.  + m ..p.p. 

*r  r rr  r r r]  ] f nyi/j  nj  1 ] 

+ ft  d jp  + F cosn.t}  (no  sum  on  r)  ( 

r vjr  r r f 

"t  h 

is  the  general  form  of  the  r normal  mode  equation  with  parametric, 
quadratic  nonlinear  inertial,  damping  and  direct  forcing  terms  included, 
ft^  - is  the  response  frequency  of  the  mode  [3]. 

Application  of  Struble's  method,  assuming  solutions  of  the  form 


pr  = + <f>r(t))  + eX± 


imposing  the  resonance  conditions  ft^  = ^2  + ^4  an<^  ^1  + ^3  = ^4  anc^  limiting 
the  analysis  to  one  of  first  order  in  e results  in  the  equations 


A3A4 

e Li  sin(K  ~ h ~ *i> 

(ftj2  - ti^2)  A3A4 

' e 4A,ft, 


' e 4A^  COS(<(,4  ' *3  - *1> 


E 4ft^  k24  Sin%  + V 


V -0)22)  A4 

~W2 + C 4A^  k24  C0S(<(,2  + V 


Al\ 

e 4 nr  l3  S1U%  - ^3 " V 

o 

( fi32“  U,32)  A1A4 

2ft^  " £ 4A^T  L3  cog(<I)4  ^3  “ 4’l) 
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k42  sin(<^2 + V - Ii4  sin(4>4 ‘ *3 " 6l)} 


(%2-^2)  ' A2  aaa3 

+ k„o  cos((J>.  + <!>.,)- 


2tt„ 


4A4«4  42  ' '^2  y4'  4AJ24 


cos(<J>4  - (j>3  - 4>x ) > 


where 


£143^32  + A134fl4*  ” (M143  + M134)i23fi4 


J'341^i2  + 2'314^42  ” ^M341  + M 314^1^4 


^413^3  + £431°1  ” (M413  + M431)niS 


and  damping  terms  have  been  excluded  for  the  time  being. 

As  can  be  seen,  modes  1 and  3 have  nonlinear  terms,  mode  2 has 
parametric  terms  and  mode  4 has  both  parametric  and  nonlinear  terms.  No 
steady  state  solution  is  expected  from  these  equations,  but  the  experimental 
evidence  points  to  'quasi-steady  state'  solution. 

Hence  as  a first  approximation  a solution  of  the  form 


4»r  = o 


A = C + C , cos(gt  + Y ) was  substituted  into  the  variational  equations 
r ro  rl  'r  ^ 


in  the  hope  of  finding  analytical  expressions  for  the  beat  frequency  3 
and  the  amplitude  ratios.  However  only  trivial  solutions  resulted,  and 
as  no  mathematical  model  was  readily  available  to  give  physical  parameters 
for  a simultaticn  of  the  variational  equations,  attention  was  focussed  on 
a simpler  two  mode  model  exhibiting  both  parametric  and  nonlinear  behaviour. 
It  was  felt  that  the  solution  of  this  simpler  system  would  give  insight 
into  appropriate  solution  forms  for  the  more  complex  case. 


4.  TWO  MODE  SYSTEM  UNDER  DIRECT  AND  PARAMETRIC  EXCITATION 


ITU  INTERNAL 


RESONANCE 


Experimental  Work 


The  photograph  on  page  13  is  of  the  two  mode  model  under  investigation 


It  differs  from  the  original  two  mode  model  in  that  bending  deflections 


occur  all  in  the  same  plane.  By  careful  adjustment  of  masses  and  lengths 


the  two  natural  frequencies  have  been  arranged  to  be  in  a 2 - 1 ratio 


Fig.  3 shows  the  natural  frequencies  and  mode  shapes  obtained 


Mode  2 


Mode  1 


Strain  gauge  pairs  at  the  roots  of  each  beam  provide  signals  representative 


of  the  generalised  coordinates  to  be  used  in  the  theoretical  analysis  of  the 


structure.  Transformation  to  normal  coordinates  is  achieved  by  using  the 


same  analogue  devices  as  in  the  four  mode  case 


The  photographs  on  pages  15,  16,  17,  18  are  of  ultraviolet  records 


taken  from  the  experimental  apparatus 


The  first,  pagel5  , shows  direct  forcing  of  the  first  mode  at  8.75  Hz 


Growth  of  the  second  mode,  through  the  internal  resonance  condition  is 


rapid 
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The  second,  page  ig,  shows  direct  forcing  on  the  second  mode  at 
17.5  Hz.  This  again  leads  to  growth  in  the  first  mode  and  a corresponding 
decrease  in  amplitude  in  the  second  mode. 

The  third,  page  17,  shows  forcing  at  the  combination  frequency,  26.4  Hz. 

Excitation  level  is  high  in  order  to  achieve  large  amplitudes  for 
nonlinear  effects  to  become  dominant  and  the  record  starts  with  the 
ordinary  forced  response,  at  the  forcing  frequency,  of  each  mode.  There  is 
then  a transformation  in  the  response  as  each  mode  is  also  excited 
parametrically  and  nonlinear  interaction  takes  place.  Work  done 
previously  by  Asmis  2 Tso  [4] , suggests  that  for  such  a system  there  is 
a continuous  interchange  of  energy  between  the  participating  modes.  This 
however  is  not  apparent  here  and  steady  state  conditions  are  achieved. 

It  should  be  noted  that  the  forcing  frequency  not  only  equals  the  sum  of 
the  two  natural  frequencies  but  also  is  three  times  the  first  natural 
frequency,  a consequence  of  the  internal  resonance  condition. 

The  fourth  photograph,  page  is,  shows  excitation  at  twice  the  second 
natural  frequency,  35.2  Hz.  This  leads  to  excitation  of  the  high  mode 
parametrically,  leading  to  growth  in  the  low  mode  through  internal  resonance 
and  a continuous  interchange  of  energy  between  the  two  modes. 

Again,  as  a consequence  of  the  internal  resonance,  the  forcing  frequency 
is  four  times  the  lower  natural  frequency  as  well  as  being  twice  the  higher. 


timing  marks  at  1 second  intervals 


second  intervals 


ks  at  I second  intervals 


timing  marks  at  1 second  intervals 
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By  approximating  the  structure  with  a two  mode  lumped  macs  model, 
with  coordinates  as  in  Figure  14 , and  applying  Lagrange's  procedure  to 
it  results  in  the  equations  of  motion: - 


{M 


q + + 2.25  y ^ )}wL  + (H^l.  5^  )w2  + - “{M^+Mgd  +2 • 2 ty 2 wq 

+ T^W2^Ws  " M2*0’9  IT^wlwl  + 0,45  £2^',v12+  1^(w2w2  + w2?)  + "l^v'lw2 

3 • 0 / * * v\ 

+ ^(W1W2  + W1W2)} 


...  1.5 

v,  )w  + 


. 1-2 
W,W,  + — S-W„W 


£1  1 1 


l2  2i 


£, 

M21,5£^l  + M2W2  + A2W2  = “ H2*(T^W2  + ~T^‘l 

1.2,  ••  . • 

- -^■(w1w2  + Wlw2)} 


where  quadratic  nonlinearities  only  have  been  retained. 

The  linearised  equations  give  expressions  for  the  natural  frequencies 
and  mode  shapes  which  can  be  adjusted  to  give  the  internal  resonance 
condition  exactly,  as  was  the  procedure  with  the  three  mode  model.  The 
equations  of  motion  in  normal  mode  form  can  then  be  formulated,  and  their 
parameters  applied  to  the  more  general  theory  which  follows. 

Equation  (1)  on  page  10  is  the  general  form  of  the  equations  of  motion. 
Denoting  the  excitation  frequency  by  the  cases  of  immediate  interest  are 

ftf  = = ^2>  + f where  - 2ft  ^ 


A solution  of  the  form  p (t)  = A (t)cos(ft  t + <j)(t))  + f cosfi^t  + ea 

*r  r r rr  r f r 

Fr 

is  assumed  where  fp  = t'  is  the  component  of  forced  response  [5]. 

r 

Application  of  Struble's  technique  yields  the  variational  equations.  The 
remainder  of  this  section  is  a catalogue  of  the  resulting  equations  for 
each  external  frequency  condition. 


Sit  • Si. 


- 2 | I A*  <°s(*z  4)  + M.  A,AacosC4- 

+ 2 d £ s.  n 4 + 2.  F,  COS  /(| 

“2/1,11,  * - 1 1 ^ f\x  ^ - M, A,Aj  sin C4l* zjO 

-2  SI*  d)((  A,  + f,  cos  rt)  + 

. 

" 2 AXA  & r C^i'  W*)  A2  " \ ?x  A,  <#  ~ 4)  + N-2  A.  Cc,sC24‘^ 


+ NYj  cx>s  <^2 


£ - ft,  ~ ^2.  fl*  +■  fA^  6lfl^ 

^ ^ A*  fl2  j 


io  here 


- f.CV1! + V2-*  - ^(^c^-AV^-xz-n-^Y 

k*  " 2 ft  A (Yu+  ' f*  A CA2  + Azi  + %+  ^zz.)  • 

A-V'v-^.x 0 

(■~^2\\  ~ mJ|f)  "^i  • 


f.  **»  ('  Ar  mzi)  ♦ Axv~  *„t)  + f.  fa  6 Z2*z‘  A»,  " "a.t  ' **».} 
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ilc  = SI, 


- £af- - ||  ^ A,  u>S2^  + n,  a,az  cos£&-**o| 


-2A,iI  = -ffp.fl,  - N,  A^S^-2^,)  -2SL^ft,: 

' * 

V 

2 fl2512  (^_  = (n’-“‘'XA;  ' fi“s&)  " ||  Nzf'X«^r'4‘)  ’ 2(2  205  4 j 
2 f\1nl  = -“£)  fi 5i,>  * 2^s;,'4  ‘2a».,lufij"j 


-2AJa1^  = 


p,  • K - * f, 

N,  = ■Si,ax^m«Li  + ,v'i.a')  “ (^>2,^1  + ^//z-^-j) 


A/j  s - ("/jii  + ^2//)-^! 


-2fi.il,  = £nL*-«D*)ft,  - | P^jMs^,  Hfi.^cps^-i^j 

-2fl.il,  - - |f  ~ N,  A, A2  C^2-  2.^,)  -2-n^A.J 

- 2A.il  '4  * (&-\-  W*)  V £\  P*fi,  C0S^,4  ^ + N^COS^-  24)] 

2 J 

- 2f)znt  - | ^ fl,  O - HZR?  sin(4-24)  -2n^d.„Azj 

P.  = ^IZ  — ^2  (J- m f*  * 4m.  fi)  ~ ^ &*,  f<  f 4a  fz)  (■•(rn..2+  mu)*  2 MU!  fi) 

?z  ■=  k2,  - il^feif,  4 44)  “ ft  12 fa  * Aif.)  + 


« 


It  should  be  noted  that  inclusion  of  the  particular  integral  in 
the  solution  leads  to  terras,  of  nonlinear  origin,  accumulating  with  the 
terms  of  parametric  origin  in  the  variational  equations. 


5.  FURTHER  WORK 

The  equations  listed  in  the  previous  section  do  not  lend  themselves 

to  analytic  solution,  and  hence  recourse  is  being  made  to  numerical 

solution  using  the  mathematical  model  based  on  the  two  mode  experimental 

model.  Some  consideration  will  be  given  to  the  effects  of  nonlinear 
\ 

elasticity,  which  are  generally  taken  to  be  small  in  comparison  with 
inertial  nonlinear  effects.  The  relevance  of  other  nonlinear  interaction 
effects  which  are  little  known  but  which  have  been  the  subject  of  research 
in  this  Department  will  also  be  considered  in  relation  to  the  class  of 
problem  under  investigation. 
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eventual  participation  of  all  the  modes,  is  presented  along  with  a theoretical 
approach  of  a more  general  nature.  Attention  is  then  turned  to  a less  complex 
two-mode  system  exhibiting  both  parametric  and  nonlinear  behaviour. 
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